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We calculate the precession of Keplerian orbits under the influence of arbitrary central-force 
perturbations. Our result is in the form of a one-dimensional integral that is straightforward to 
evaluate numerically. We demonstrate the effectiveness of our formula for the case of the Yukawa 
potential. We obtain analytic results for potentials of the form V(r) = or" and V{r) = aln(r/A) 
in terms of the hypergeometric function 2F1 (5 — -j, 1 — -j; 2; e^), where e is the eccentricity. Our 
results reproduce the known general relativistic (n = —3), constant force {n — 1), and cosmological 
constant (n = 2) precession formulas. Planetary precessions are often used to constrain the sizes of 
hypothetical new weak forces-our results allow for more precise, and often stronger, constraints on 
such proposed new forces. 

^ ' PACS numbers: 04.90.-|-e, 98.80.Jk 
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. By the early years of the twentieth century, several workers had calculated the expected precession of the orbit 
of the planet Mercury based on the equations of Newtonian mechanics taking into account the perturbing effects of 
other planets in the solar system. (See Doolittle [l| for a discussion of methods, a summary of results, and early 
J> ■ references, and Stewart [2| for a recent discussion.) The prediction was for a precession of 575 arcsec per century (in 
^~f^ ' a sun-centered coordinate system). Observation found the actual precession to be 532 arcsec per century, leaving 43 

■ arcsec per century unexplained. This discrepancy remained a mystery until 1915, when Einstein, in the course of the 
I development of general relativity, calculated the weak-field expansion of his new theory and predicted an additional 

■ precession of just this amount. ^ Since then, the precession of Mercury and the other planets has been used 
— . ' to place bounds on the size of additional hypothetical forces. Examples of such forces include the "anti-gravity" 

. due to the cosmological constant [1, [1, 0] , forces associated with a distance dependence of the gravitational constant 
Q , forces proposed as alternatives to dark matter [l3| , forces proposed as phenomenological explanations of the 
Pioneer Anomaly [ll^ 12j_ 13^3, [l^, and forces induced from higher dimensional models [16, 17, 18, 19] or extensions 



I. INTRODUCTION 



^ ' of general relativity 



^jrj' In most of the precession calculations, the assumption of nearly-circular orbital motion has been made in order 
• • to simplify the calculations. We show in this work that such an approximation is not necessary, and give a simple 
. !^ formula for the calculation of the precession induced by any central force even for highly eccentric orbits. Our result 
, can be used to get numerical values for the precession in all cases, and analytic values for power-law and logarithmic 
^ ' potentials. Since the exact precession values can differ from the nearly-circular values by a factor of two or more, 
5^ , significantly more precise results for constraints on proposed new forces can be obtained through use of the present 
work. 

The outline of this work is as follows. In Sec. II we review the case of nearly-circular orbits and give a derivation of 
the relevant precession formula. In Sec. Ill we develop a formula for the precession of non-circular orbits in the form 
of a one-dimensional integral. In Sec. IV we apply our formula to obtain precession results for a perturbing force of 
Yukawa form as a function of the Yukawa range parameter and the orbital eccentricity. In Sec. V we obtain analytic 
results for power-law and logarithmic potentials. Finally, in Sec. VI we discuss a number of examples of hypothetical 
central-force perturbations for which precession bounds are appropriate. 

II. PRECESSION IN NEARLY-CIRCULAR ORBITS 

We begin with a brief review of the Keplerian case. The unperturbed orbit equation is 

(fu GM 
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where r = \/u and </? are cylindrical coordinates in the orbital plane, G is the gravitational constant, M is the central 
mass, and h — ip/u^ is the angular momentum per unit mass of the orbiting body. The standard elliptical solution 
has the form 

^*o(<^) = {1 + ecos((^ - (/3o)} (2) 

where L = h'^/{GM) is the semilatus rectum of the orbital ellipse and e is its eccentricity. The maximum and 
minimum distances from the orbiting body to the force center have the values 

^± = T^- (3) 
The semimajor axis a is thus related to the semilatus rectum according to 



«=^(^++^-) = (4) 



so that r± — {IzL e)a. The period r is given by Kepler's Third Law: t'^ — Att'^o^ /{GM). 

In the limit of nearly circular orbits we can find an approximate solution to the orbital equation that displays the 
precession of the axes of the ellipse. In this approach we generalize the method of Ohanian [23J. We assume the 
presence of a radial perturbing force F{r) in addition to the usual Newtonian gravitational force —GMm/r'^. The 
perturbed orbit equation takes the form 



(Pu GM giu) 
\- u — — 



,2 ■ " = J^T- (5) 



where g{u) = r^^^-^\r=i/u- The unperturbed solution is given in ^ above. We assume that the eccentricity e is 
small, and expand 

GM 

giu) = ga+ gi-^ecos{(p ~ tpo) + 0{e'^) (6) 

where go — g{l/L) and gi — g'{\/L). To first order in e the orbit equation can be written as 

dFu GM - ffo GM 

^+"~ 5i-^ecos(^-(^o). (7) 

We try a solution 

u{^) = { ^^\^ {1 + ecos [(1 + B)[^- (/.o)]} (8) 

and find, to first order in the perturbation, that B = gi/{2h'^). It follows that the precession angle (per revolution) is 
A0p = —27:3 = —ngi/h^. We notice immediately that a perturbation, having g{u) = const., gives no precession. 
If we write the perturbing force in terms of a potential 

Fir)^-^, (9) 
ar 

we find that 

d ( 1 2dV\ . 1 d^V. 

au \ m dr J m du 

where we have expressed the potential V(r) in terms of u = 1/r. It follows that the precession /S.9p — —ngi/Ii^ can 
be written in terms of the potential as 

71- Id'^V, , , 

^'^--GM^-L^\--^I^ 

where we have used Y? = GML. This is equivalent to a formula for the nearly-circular precession that has been used 
by Dvali, Gruzinov, and Zaldarriaga [l7j . 
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For a power-law perturbing potential V{r) = anr"' the results are particularly simple. The force has the form 
F{r) = — na^r""""^, and the precession (per revolution) is 

Ae» = -n(n + l)7r^^. (12) 

There is an interesting relationship between the precession angles for powers n and — (n + 1). The precession when 
n — > — (n + 1) is 

A9,{-{n + 1)) = -n(n + l)n , (13) 



so that 



A^p(n) ^ ^2n+i Agp(-(n+l)) _ 



(14) 



III. PRECESSION IN NON-CIRCULAR ORBITS 



In this section we will obtain an expression for the precession valid for orbits of any eccentricity and explore some 
of its properties. In the following sections we will discuss the Yukawa force as an example where numerical techniques 
arc required, and wc will find analytic results for a number of useful potentials. 

We will use energy conservation to obtain a simple expression for the precession. The conserved energy for motion 
in potential —GMm/r + V{r) has the form 



m 



GMm 



E=^[f-' + '-^]-^ + yir). (15) 



In terms of the more useful variable u = l/r this is 



E=^\^\^—j+u']-GMmu + V{u). (16) 



We solve for du/dcp and integrate to obtain 



{ '2E I 2GMu _ .,2 _ 



where the sign is determined by the position in the orbit (increasing radius or decreasing radius). 

Returning to the non-pcrturbed problem for a moment, we see that the turning points are the solutions a± = 
1/rip = (1 ± e)/L to the equation 

2E 2GMu 2 / X/ N 

^^^^~h^~'' ^^^^ 

The angle change for a full orbit (twice the change for the interval u = a- to u = a+) is 

r"* du 

^/{a+-u){u-a-) 



r 

Aip = 2 1^ , = 27r, (19) 



as expected. Explicit expression for o± are not needed in the evaluation of this integral. 

Wc arc now prepared to write down the general expression for the angle increase per complete radial oscillation (u 
from smallest to largest and back to smallest). It is 

= 2 r —TT. ' (20) 



«- ( 2-E I IGMu „.2 2V(m) ^ 



where u± are the zeros of 



2E 2GMu 2 2F(u) 
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with u± « a±. The precession per revolution is thus A9p = A(p — 2tt. 

We will need approximate values for the zeros u±. We note that the expression D{u) can be written as 

D{u) = {a+-u){u~a_)~'^:^^. (22) 
We write u± — a± + 5±, and find that the first-order shifts in the positions of the zeros are given by 

°± = 7 ■ (23) 

It is useful to express D{u) with its zeros explicitly displayed as 

D{u)^{u+-u){u^u^)(^ + ^^ . (24) 

It is clear that Q is small, of the order of the perturbation, since D{u) can also be written as in (I22p and u± « a±. 
An explicit approximation for Q can be found by equating the two forms for D{u). We find that 



-(a+ - a-)V{u) + (a+ - u)V{a-) + {u- a-)V{a+) 
(a+ — a_)(a+ — u)(u — a_) 



(25) 



We note that the numerator of Q vanishes when a_|_ = a_ and when u — no Q \& non-singular for the whole physical 
range < u < a_|_. With this expression for D{u), the precession can be written as 



du f 2Q 



-1/2 



(26) 



\/(a+ - - a-) 



A number of consistency checks can be applied to this result. First, we notice that Q vanishes when either V{u) — 
const, (no force) or V{u) cx u ( a contribution to the force), so there is no precession in these cases. Also, we can 
work out the precession for a nearly-circular orbit as a limiting case of (|26|) for a+ — > a_. When a+, a-, and u are all 
similar, we can use a series expansion to show that Q {1/2)V" (u). The corresponding precession, from (|26p . agrees 
with our earlier result (|lip . For a power-law perturbing potential V{r) — a„r", the change of variable u = a^a-/x 
relates A6'p(n) and Mp{-{n + 1)). We find that 

A0p(n) _,2n+iA0p(-(n + l)) ^27) 



an a- 



(«+i) 



which is the generalization of (jl4p to non-circular orbits. Here 6 = a\/l — = L/Vl — = is the 

semiminor axis. 

Occasionally it is more convenient to express the perturbation as a force instead of a potential, for instance when 
the phenomenological perturbing force has a complicated structure or is given as an anomalous acceleration. Our 
precession result (1^5)) can be rewritten in terms of the force by means of an integration by parts (using J du{{aj^ — 
u){u - a_))"^/^ = 2(2u - a+ - a_)(a+ - a_)~^((a+ - u)(u - a_))^^/^) in the form 

M = -4 r du{2u-a+-a^)V'{u) 

^ GMmL{a+- a^f ^ {a+ - u){u - a^) 

An additional advantage of this form is the absence of the potentially troublesome factor ((a+ — u)(u — a_))^/^ in 
the denominator. For practical use it is convenient to use a dimensionless integration variable with a fixed range. We 
define m = (1 + ez)/ L, and find 

_ -2L^ dzz F{z)lm 

where F(z) is the perturbing force at radius r = L/(l + ez). An alternative expression for the precession, in terms of 
the potential, is 

AO C _Jl^dViz) 
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where V{z) is the perturbing potential at radius r = L/{1 + ez). 

A general consequence of either form or ([50]) is the fact that the precession AOp contains only even powers of 
the eccentricity e. This is true because only even powers of z contribute in the integrals, and e and z are tied together 
in the definition of r. 

Of course, celestial dynamics is an exceedingly well-developed discipline, and expressions equivalent to ((29|) and 
()30|) can be derived from the formulas of orbital perturbation theory given in the standard texts (see for example 
Danby Landau and Lifshitz [25*1 give a precession result in an alternative, but also equivalent, form. We have 

found the expressions (^9]) and (PUI useful for further developments. 



IV. THE YUKAWA FORCE 



— r/X 

V{r) = a- (31) 



The Yukawa potential 



of strength a and range A is a field-theory-motivated [261 127| assumption for a correction to the gravitational force. 
In this section, we work out the precession due to a perturbation of Yukawa form. 

The precession A0p(K, e) due to a Yukawa perturbation depends on two parameters: a range parameter k = L/X 
and the eccentricity. It seems difficult to obtain an analytic form for A0p(K, e), so we work numerically instead. Use 
of either or ([50]) leads to the integral form 

The small-e limit of this, obtained either as the first non- vanishing term in an expansion in e or by use of (|lip directly, 
is 

A^p(«^,0) = ^^^2e-''. (33) 

We take A9p{K,0) to set the scale for the Yukawa precession, and focus instead on the effect of eccentricity encoded 
in the relative precession I{k, e) where 

A6ip(K,e) = A6'p(K,0)/(K,e). (34) 

The integral for the relative precession is 

2 dz z / n \ ( nez \ 

Iin,e) = ^ l + exp . 35 

TTK-^e VI - ^2 V 1 + ez/ \l + ezj 

The relative precession can give an important correction to the zero-eccentricity result. For example, the relative 
precession of the asteroid Icarus (with eccentricity e=0.827) with a Yukawa perturbation having range parameter 
K = 0.1 would be 1(0.1,0.827) = 4.57. The general behavior of I{K,e) as a function of its parameters is shown in 
Fig. 1. For most values of its parameters I{k, e) > 1, and the constraints on the Yukawa strength parameter a from 
precession bounds are strengthened by taking the eccentricity dependence of /(k, e) into account. 



V. ANALYTIC RESULTS FOR POWER-LAW AND LOGARITHMIC POTENTIALS 

For some perturbing potentials it is possible to obtain analytic precession results. In this section we work out the 
analytic results for power law and logarithmic perturbations. There are known expressions for the precession due to 
a few particular power law perturbations. The small-field expansion of general relativity, valid outside a spherically 
symmetric mass, leads to a potential proportional to r~^. A constant force has a potential proportional to r. Also, a 
uniform dark energy force is modeled by a potential proportional to . We give a general expression for a perturbative 
potential V{r) = anr^ that contains these results as special cases. 

For a power law potential we take V(r) = anr"^ and apply (|30p . It turns out to be easiest to work first with a 
negative exponent. In particular, for a potential V{r) = a_(„+i)r^("+^' — a_(„+i)(l + ez)"+^/L"+^, the precession 
is 



-2a_(„_|_i)(n -I- 1) dzz 

z'- 
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FIG. 1: Contour plot showing the relative precession I{n,e) for the Yukawa potential as a function of the range parameter k 
and eccentricity e. The contour lines are at 0.75, 0.85, . . . , 1.55. The darkest (lightest) shading corresponds to the smallest 
(largest) values for the relative precession. The medium gray shown on the left side of the plot corresponds to values near 
one. As expected, the largest variations from one occur for large eccentricity. Surprisingly, there is a range of k for which the 
relative precession is depressed from its e = value. 



The integral can be done by expanding (1 + ez)" in a series and integrating term by term. One finds that 

oo 



1=0 



i=0 



2i + V (2z + 2)!! 



In the first step here we made use of the fact that only even powers of z survive the integration. Our complete result 
then for negative powers is 

where 

X„(e) = „(n + 1) 2^1 Q - |, 1 - ^ ; 2 ; e^^ . (39) 

The precession for positive powers can be obtained by using (|27p or the hypergeometric identity (28l | 

X„(e) = (l-e2)"+V2^_(„+i)(e), (40) 

and is 
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TABLE L Polynomials Xn{e) for various values of n. 



n Xn(e) 

1 2 

2 6 

3 3(4 + 6^) 

4 5(4 + 36^) 

5 15(8 + I2e2 +e*) 

6 %(8 + 20e2 + 5e'*) 



7 1(64 + 2406^ + 120e'' + 5e®) 



We have chosen to write A9p in terms of the semilatus rectum L for negative powers and in terms of the semimajor 
axis a for positive powers in order to simphfy the formulas. We note that expressions ([55]) and (|4ip are in fact both 
valid for all n. 

The hypergeometric function 2F1 terminates whenever one of its first two arguments is zero or a negative integer, 
which happens for any natural number n. More generally, 2^1 and A9p(n) are well defined for any n, integral or not. 
Results for the first few polynomials Xn(c) are given in Table I. 

Our power law expressions can be put to immediate use to reproduce known results for n — —3, n = 1, and n = 2. 
Einstein Q considered orbits around a spherically symmetric mass in the weak- field limit. The first post-Newtonian 
correction is a perturbing potential V{r) = — ^^^2^^ where h? = GML. The corresponding precession from ([55]) is 

, , , GnGM , , 

A0p(GR) = . (42) 

This gives the famous 43 arcsec per century when applied to the orbit of Mercury. A constant perturbing radial force 
F = mao is represented by the potential V{r) — —maor. The associated precession (per revolution) from (|4ip is 

27raoa^ 



AOpiF = mao) = ^1^^ (43) 

consistent with the result of Sanders [22] . The extra-Newtonian radial force due to a uniform dark energy is given by 

M 

F{r) = -mAc^r . (44) 
3 

This corresponds to a perturbing potential V{r) = — ^mAc^r^. The corresponding dark energy induced precession 
(per revolution) from (|4ip is 

A0p(A) = ^^v/r^ (45) 

in agreement with the expression of Kerr, Hauck, and Mashhoon [30]. (We note the curious fact that the n — —3 
general relativity precession and the n = 2 dark energy precession are dual to each other in the sense of 
(|27)) .) Result (|45|) can be used to find a bound on the cosmological constant A. The measured anomalous precession 
for Mercury, after accounting for known Newtonian and relativistic effects and solar oblateness, is —0.0036(50) arcsec 
per century [si^. The orbital parameters for Mercury are a = 5.79 x 10^°m, e = 0.206, r ~ 0.241y = 7.60 x lO^s, and 
the solar mass is M = Mq = 1.99 x 10'^'^ kg, leading to a bound on A of 

- 2.5 X lO-'^'^cm"^ < A < 0.4 X 10~'^^cm-'^ . (46) 

Tighter bounds on A can be found by using precession data for other planets. ^ These bounds are several orders of 
magnitude larger than the currently accepted value of A found from cosmological and large-scale measurements, but 
it is interesting to find bounds on the cosmological constant based on purely solar system measurements Q . 

An analytic result can also be obtained for a logarithmic potential. We assume a perturbative potential modification 
of the form 



V{r) ^ aln(J^^ (47) 



8 



where A is an arbitrary scale parameter. Then V'{z) = —ae/{l + ez), and from ((30l) the associated precession is 
We have already done the required integral in ([57)1 : 

1 dz Z TT „ 3 „ 2^ 



where [33 | was used to obtain the value of the hypergeometric function. The precession for a logarithmic perturbation 
is 

We note that A0p(log) = hm„^o 



VI. SUMMARY AND DISCUSSION 



We have found expressions for the precession induced by central-force perturbations that generalize the known 
nearly-circular result (fTT|) . The most general results are ((^ for perturbing force and ([50]) for perturbing poten- 

tial y^z\ These expressions are straightforward to evaluate numerically. For power-law and logarithmic perturbing 
potentials, simple analytic results can be found: ([55)1 and (j3T|) for power-law potentials, and (|50p for a logarithmic 
potential. As examples of the use of these formulas, we have shown in some detail how the precession for the Yukawa 
potential can be found numerically. Also, we applied our formulas to reproduce the known precession results for 
the general relativistic n = —3 potential, the constant force n = 1 potential, and the cosmological constant n = 2 
potential. 

Many other central-force modifications to gravity can be found in the literature, all of which can be treated according 
to the methods developed here. A small sample of proposed modifications to gravity that have been discussed in the 
recent literature include the following. 

There are many proposals for modifications to gravity in our four-dimensional world that arise from a more funda- 
mental theory in a higher dimension. An example is the five-dimensional braneworld scenario of Dvali, Gabadadze, 
and Porrati [l6[, in which the short distance limit of four-dimensional gravity has the form 



V{r) = 



-GMt 



2 

1-f - 

TT 



To 



7-1 



-)+0(r2) 



(51) 



Further study of this model by Gruzinov [Toj and by Dvali, Gruzinov, and Zaldarriaga |T?| points to a perturbing 
potential proportional to r^^^. Capozziello et al. [23| discuss the effect of general Yukawa corrections 



V{t) 



GMm 



E 



-r/Xk 



(52) 



Capozziello and Lambiase [l8[ studied the Newtonian limit of string-dilation gravity, and found a perturbing potential 
containing a linear combination of Yukawa, r^, and cosh(r/A) terms. 

The Pioneer Anomaly, the anomalous acceleration of about 8.5 x lQ~^^m/ for 2QAU < r < 70AU directed towards 
the sun experienced by the Pioneer 10/11 spacecraft, has been a strong source of motivation for the development of 
modified theories of gravity. One model of the proposed gravitational force considered by Anderson et al. [ll|, [l2| is 
a Yukawa potential (Newtonian plus perturbation) that has the usual Newtonian 1 /r^ force but no 1 /r correction: 



V{r) = - 



GMm 1 



1 



(53) 



Jaekel and Reynaud [l3| consider models involving linear and quadratic perturbing potentials. Brownstein and 
Moffat [l3| discuss the Newtonian limit of the scalar-tensor vector gravity theory [2l|, which is best expressed as an 
acceleration of Yukawa-like form but with scale-dependent parameters: 



F{r) 



m 



GM 



1 + a{r) 



1 



A(r) 



(54) 
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For any explicitly given strength and range parameters a{r) and A(r), the precession due to such an anomalous 
acceleration can be evaluated numerically using (P^j) . Sanders [l^ discussed solar-system constraints on multi- field 
extensions to general relativity that might explain the Pioneer Anomaly. 

For order-of-magnitude estimates of orbital precession it is often adequate to use the near-circular formula (jlip . 
In some situations, especially for highly eccentric orbits, and whenever precision is desired, precessions should be 
calculated taking into account the eccentricity of the orbit. This is straightforward to do numerically for any perturbing 
potential or force, and for power-law and logarithmic perturbing potentials analytic results are now available. 
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